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Abstract. Global degree/strength based preferential attachment is widely used as 
an evolution mechanism of networks. But it is hard to believe that any individual can 
get global information and shape the network architecture based on it. In this paper, 
it is found that the global preferential attachment emerges from the local interaction 
models, including distance-dependent preferential attachment (DDPA) evolving model 
of weighted networks(M. Li et al, New Journal of Physics 8 (2006) 72), acquaintance 
network model(J. Davidsen et al, Phys. Rev. Lett. 88 (2002) 128701) and connecting 
nearest-neighbor(CNN) model(A Vazquez, Phys. Rev. E 67 (2003) 056104). For 
DDPA model and CNN model, the attachment rate depends linearly on the degree 
or vertex strength, while for acquaintance network model, the dependence follows a 
sublinear power law. It implies that for the evolution of social networks, local contact 
could be more fundamental than the presumed global preferential attachment. 
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1. Introduction 

Power law distribution commonly exists in complex networks, including planned network 
such as Internet p. j and unplanned social networks, such as scientific collaboration 
networkspj El HJ [5], actor collaboration networks [§], peer-to-peer networks [Z], mobile 
networkspB] and much more. Their degree distributions obey power law, and if 
applicable, vertex strength(summing the weights of links that connect to a vertex) [U [5] 
and link weight [9 J also follow power law. Many models have been proposed to capture 
the topological evolution of complex networks (see reviews [101 [II]). Especially a 
class of models based on the idea of global preferential attachment, firstly proposed 
in BA model [12], is quite successful to reproduce the power law distributions of 
degree/strength [9l [12j [131 E3]- From the degree-based preferential attachment in BA 
model, later Barratp2] developed a model of weighted networks with strength-based 
preferential attachment. 

Some networks have hubs or data centers, where global information is collected for 
future use. For example, WWW has search engines and Internet has routers and DNS 
servers. For networks with such centers, global preferential attachment mechanism 
may seem reasonable. However, in social network there is no data centers collecting 
and providing global information. It is impossible for every individual to know global 
information of system, such as degree or strength of every individual. For social 
networks, evolution mechanism based on local quantities would be more natural. 

Actually, there are already several evolving models for scale-free networks based 
on local rules, such as distance-dependent preferential attachment (DDPA) model[15j, 
acquaintance network modelJH], connecting nearest-neighbor(CNN) model [IT], random 
walk modelpSl HE], redirection modeljSHl EI], optimization model[22j [23] and so on. 
For instance, in DDPA model a measure of closeness relation is defined locally within 
second neighbors and then links are built up preferentially according to the relation. 
Acquaintance network model evolves via people introduced to know each other by a 
common acquaintance [H]. In CNN modelflT]. at every time step, a new vertex is added 
or a potential edge within second neighbors is converted into an edge. All above models 
based on local rules can also reproduce the common topological characters of complex 
networks. 

We notice that both the two potentially conflicting types of models successfully 
capture the main characters of scale-free networks, but one needs global information 
while the other only requires local information. Here we suggest a wild guess to resolve 
the conflict and make the whole picture more consistent. We conjecture that global 
preferential attachment emerges from local contact based models. In this paper, we are 
going to investigate whether or not the above conjecture holds in the following local 
models: DDPA modelpj)], acquaintance networks model[l6] and CNN model[TT]. 

This paper is organized as follows. In section EJ we first briefly describe the method 
used to check preferential attachment. In section [3], we apply this method to DDPA 
model, acquaintance networks model and CNN model, and present the results. At last, 
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in Section H] we give some concluding remarks. And it turns out our conjecture does 
hold in the above models. 

2. Methods for Measuring Preferential Attachment 

The preferential attachment hypothesis states that the rate H(k) with which a vertex 
with k links acquires new links is a monotonically increasing function of /c[12], namely 



n(*,) = =^ = cm. (i) 



For BA model a = 1[12]. 

H. Jeong et al[2I] proposed a method to check global preferential attachment from 
data on evolution process. Here we will briefly describe the method. In evolving process, 
we record the order of each node and link joining the system within a relatively short 
time frame after the network evolves steadily in a long time. At a large enough time 
To, consider all vertices existing in the system, called To vertices. Next select a time Ti 
(Ti > T and (T x — T ) <C T ), add a new vertex at every time step between [T ,Ti]. 
Firstly, count the number N(k) of such vertices with exactly k degree in T vertex. 
Secondly, record all the vertices in T vertex to whom the new links are attached as 
Q. Lastly, count the number of vertices with exactly k degree in Q as A(k). At this 
condition, Il(/c, T , Ti) will be independent of T and Ti but depend on k only [24]. 

A convenient definition of H(k, To, Ti) function could be defined as, 

#) 

n(*,r 0l r 1 ) = -^. (2) 

^k' N(k') 

According to mean field theory and Eq.flTJ, A(k) = M^f^i, k)kf / J2j kf), where M 
is the number of new links attaching to vertices Q. So Eq.Q can be written as 

n(fc) = ^4 = U = (3) 

EAME, S(kj, k>)^/N{k>)) £* £* k'° 

where J2i°~(ki, k)kf = k a N{k) is used. To avoid the effect of noise, we also study the 
cumulative function instead of H(k) 

n(k) = f k Tl(k')dk'. (4) 
Jo 

If there is a global preferential attachment, then 

K(Jfe) oc k a+1 . (5) 

In order to prove the validity of formula (J2J), we apply this method on BA model, 
starting from a fully connected Uq initial network, where (k) = 2m, a = 1 and 
7 = 3. The exponent in Eq.()3]) should be 1 and the measurement shows that they 
are approximately 1 (as shown in Fig. [1]). This indicates that the formula (EJ) is fit to 
measure the form of U(k). 
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Figure 1. The n(k) function determined numerically for BA model based on Equation 
([2]). The slope of solid line is 2 (corresponding to a = 1 as a + 1 = 2), and the slope of 
dash line is 1 (a = 0). This indicates that the exponent a is almost 1. If not mentioned 
in the following figures, we run the model for 100 times with the same parameters 
to count A(k,To) and N(k,To), and then we determine K(fc)(fc(s)) according to Eq. 
([2]) and estimate the value of a, where the power-law exponents are calculated on 
intermediate degree (vertex strength). We used T\ = Tq + 100 if not mentioned. The 
parameters are no = 10 and m = 5. 



In measurement, we only focus on the vertices new links are attached onto, so we 
only consider ending vertices. Furthermore, measurements on external links from the 
new vertex and internal links among existing vertices, can be done separately when 
necessary. 

3. Measuring preferential attachment of Local Contact Model 

In this section, we are going to check the possibility of such global attachment emerging 
from the following three local contact models. 

DDPA model: The evolution of networks starts from a fully connected no initial 
network. At every time step, one new vertex is added into the network by connecting it 
randomly to one old vertex. And then other I old vertices are randomly activated. Every 
one (denoted as vertex n) of these 1 + 1 vertices can attempt to build up m connections. 
The probability for every link (except the first link from the new vertex) from a vertex 
n to a vertex i is given by 



where d\ means neighbors of vertex n up to the second separation. L ni is the similarity 
distance [3 |9], which means that the larger is the similarity distance, the closer is the 
relation between the two end vertices. If two vertices i and j are connected directly, 
Lij = Wij, where ib^ is the similarity weight, meaning the larger the closer, such as the 
number of collaboration in scientific collaboration networks [2j [3l HJ E], and the duration 
of calls in mobile networksjH]. Otherwise, is calculated as follows. If vertex i (fi) 



n 



in 



(6) 
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Figure 2. The n(k) (a) and k(s) (b) function determined numerically for external links 
of DDPA model. The insets are results for the first external links, where = 0.01 
and a s = 0.05, almost flat distributions. And (a*) = 0.85(a) and (a s ) — 0.91(b) for 
all other external links. The parameters of DDPA model are no = 10(the size of initial 
networks), I = l(the number of activated old vertex) and m = 5 (the number of new 
connections on every activated vertex). 



is connected to vertex \i (j) with similarity link weight w ifl then L\j = 




The final similarity distance between vertex i and vertex j is = max jL^j, the 
maximum one via all vertices {fi}. In other words, it is sufficient to know the degree 
of familiarity between second neighbors. In addition, the reconnection between linked 
vertices is allowed. When building up connection between vertices i,j represents an 
event happening between the vertices, the number of occurrences of such an event may 
be defined as the connection count 7^-, which will be transferred into similarity link 
weight Wij = f(Tij), e.g., Wij = Tij. This DDPA mechanism does not use vertex degree 
or vertex strength, i.e., Sj = J2j w iji as the reference, not even locally. It only makes use 
of the information about local similarity distance within second neighbors, but scale-free 
phenomenon does appear in DDPA model for intermediate degrees/strengths(see Fig.6 
in Ref.pj)]). 

For DDPA model, we need to do this measure for three kind of links: the first 
external links from the new vertex, all other external links from the new vertex and 
internal links. The same method can also be applied to check preferential attachment 
according to vertex strength. Since the ending vertex of first external link is selected 
randomly, the probability of a vertex selected should be independent of its degree (k) 
or strength(s). In the inset of Figj2j, we see the first external link does follow flat 
distribution. But for all other external links, we find that K(k) (measuring according to 
degree) and k(s) (measuring according to strength) for intermediate degrees/strengths 
follow a straight line on a log-log plot, indicating IT(fc) cx k ak or Tl(s) oc s as . The curves 
of n(k) for DDPA model are consistent with those for BA model (FigJTJ). And they 
are parallel with different To (as shown in FigEJ), indicating that H(k) and IT(s) are 
independent of T and 7\, and depend on k or s only. Their power law exponents are 
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Figure 3. The K,(k) (a) and n(s) (b) function determined numerically for internal links 
of DDPA model. The curves increase faster than linear(shown as dashed line in the 
figures) , showing that the probability of vertex linked by internal link is proportional 
to degree and vertex strength with (ctk) = 1.07(a) and {a s ) = 1.1(b). Inserted are 
the values of (a) for different choice of T . It shows very small variation. 



respectively (at&) = 0.85 and (a s ) = 0.91. 

In the measurement of internal links, we find that global preferential attachment 
is also valid. The curves of n(k) and k(s) follow a straight line on a log-log plot. Their 
power law exponents are respectively = 1.07 and (a s ) = 1.1 (as shown in Figj3]). 
Actually, if checking all links together, including the first external link, the curves of 
K,(k) and k(s) follow a straight line on a log-log plot, where exponents are respectively 
{dkj = 0.98 and (a s ) = 1.04 when 1 = 1. The parameters almost have no effect on the 
exponents (a*.) and (a s ). 

Acquaintance model: At first, a random network with N vertices, where every 
pair vertices are connected with probability pi, is generated. Then acquaintance 
networks evolve according to the following rules: (i) One randomly chosen person 
introduces two random acquaintances of his to each another. If they have not met 
before, a new link between them is formed. If he has less than two acquaintances, then 
he introduces himself to another random person, (ii) With probability p one randomly 
chosen person is removed from the network, including all links connected to this vertex, 
and replaced by a new person with one randomly chosen acquaintance [H]. The model 
generates degree distributions spanning scale-free and exponential regimes(see Fig.l in 
Ref. [IS]). Though the size of network is fixed, this model includes birth and death 
process. New links emerge between individuals frequently. After evolving for To times, 
again we record and analyze evolution process between T and Ti. We can see that n{k) 
also follows a straight line on a log-log plot, where a ~ 0.80(as shown in Fig. H|). 

CNN model: Next we consider a variant of the acquaintance model, connecting 
nearest-neighbor model [17] with increasing network size. CNN model starts from a 
single vertex and evolves according to the following rules: (i) With probability 1 — u, a 
new vertex(denoted as n) is introduced into the network by connecting it randomly to 
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Figure 4. The n(k) function determined numerically for acquaintance network model. 
Inserted are the values of (a) for different choice of To. Results from fitting of a shows 
that (a) = 0.815. The parameters are N = 10000, pi = 0.01 and p = 0.01. 




Figure 5. The n{k) function determined numerically for all links of CNN model. 
Inserted are the values of (a) for different choice of T . For u = 0.4 we find that 
(a) w 1.00. 



one old vertex j. The pair of vertices [n, i] between all nearest neighbors i of vertex j and 
the new vertex n is recorded in d([n,i]). (ii) Then, with probability u two vertices [n,i] 
selected randomly from d are connected. When the size of network is sufficiently large, 
the distribution of intermediate degrees exhibits a power-law decay P(k) ~ /c~ 7 (see Fig. 
8 in Ref. p2]), and the value of 7 depends on the parameter u. 

In the process of measurement, we record the degree of vertex j and the other two 
selected end vertices. We see that n(k) follows a straight line on a log-log plot, and they 
are parallel with different T (as shown in Figj5]). For different parameter u, the value of 
exponent a are different, for example, a ~ 1.11 for u = 0.3 and a ~ 0.95 for u = 0.7. 
In Figj5j we show the measurement on u = 0.4. And there we can see that a ~ 1.00. 

Although a is close to 1 for DDPA and CNN model and a = 0.80 for acquaintance 
model, we find that n(k) ~ k a+1 holds in all models. And this indicates that although 
links are created by local information, it seems as if the network evolves according to 
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global preferential attachment. 
4. Concluding Remarks 

In this paper, we measured global degree/strength preferential attachment of DDPA 
model, acquaintance network model and CNN model. These models make use of 
local quantities instead of global degree/strength preferential attachment mechanism. 
However, our measurement shows that they still can be seen as if the networks 
evolving according to global preferential attachment. From this point of view, 
preferential attachment is an emergent phenomenon from the more fundamental local 
rules. Empirical study on e-mail networks in Ref. [25] also suggests that cyclic closure, 
especially triadic closure, plays an important role in building up new links. In modeling 
social networks, this observation may suggest that local rules are preferred rather than 
the global degree/strength preferential attachment. For people who believe in the later, 
our discovery is especially meaningful. Our results indicate that one does not need to 
worry about the widely believed and used global preferential attachment because it can 
emerge from properly designed local rules. 
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